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FOREWORD 


This  document  is  the  first  scientific  report  on  AFGL  Contract 
FI9628-76-C-O296 . The  work  described  herein  represents  the  completion 
of  the  following  contract  items: 

a„  Modify  the  electron  density  routine  of  the  existing  ^-D 
ray  tracing  program  to  accept  a standardized  trough  wall 
profile  along  an  arbitrary  smoo^h  curve  given  in  geographic 
coordinates  = 

b Using  the  trough  wall  model  provided  by  LIB,  compute  the 
effects  of  : to  5 trough  wall  configurations  describing 
median  conditions,  and  moderate  and  extreme  deviations 
thereof,  on  frequencies  in  the  5 to  lo  MHz  range.  A 
minimum  of  four  bores ights  for  a narrow  beam  antenna  will 
be  selected  for  the  computations.  The  range /az imuth 
errors  and  the  relative  energy  density  of  the  propagated 
beam  in  the  landing  area  will  be  determined  and  shown  as 
a function  of  boresight,  frequency,  and  chosen  trough  wall 
model . 

Use  was  made  herein  of  work  performed  under  AFCRL  Contracts  F19d28-7J" 
C-0j07  and  FlSi628-76-C-0029 . 
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Introduct ion 


I . 

The  presence  of  the  auroral  oval  with  its  sharp  horizontal 
gradients  especially  in  the  region  of  the  poleward  F-layer 
trough  wall  causes  large  azimuthal  deviations  from  the  great 
circle  path  for  radio  waves  in  the  HF  band.  It  furthermore 
drastically  changes  the  range  of  a ray  to  its  landing  point,  if 
compared  with  a case  of  a homogeneous,  horizontally  stratified 
ionosphere.  For  the  case  of  an  OTH  backscatter  situation  this 
results  in  both  range  and  azimuth  errors  of  the  detected  target 
unless  corrections  for  the  effect  of  the  gradients  are  included 
in  the  target  registration  process.  Using  the  AFGL  three- 
dimensional  ray  tracing  computer  program  (Langworthy  and  Barrett, 

1 'T'p)>  various  cross-sections  derived  from  a new  trough  wall  model 
have  been  simulated  and  backscatter  results  have  been  tabulated 
to  give  an  indication  of  anticipated  deviation  in  ground  range  and 
azimuth.  The  nower  with  which  these  signals  may  be  returned 
relative  to  signals  not  experiencing  effects  of  gradients  associated 
with  the  edge  of  the  auroral  oval  has  also  been  determined. 

The  data  used  to  simulate  the  trough  wall  throughout  this 
report  is  based  on  a model  derived  frcmAlouette  topside  ionograms . 


2 . Ionospheric  Model 

In  order  to  determine  the  magnitude  of  range  and  azimuth 
errors,  which  an  OTH  system  will  experience  if  operating  in  the 
vicinity  of  the  poleward  F-layer  trough  wall,  yD  ray  tracing  technique  was 
iipplied  to  a rew  model  of  the  trough  wall  (Pike,  1976).  The  location 
of  the  equatorward  edge  of  the  poleward  trough  wall  is  specified 
by  a series  of  points  in  the  geographic  coordinate  system.  These 
points  are  then  converted  to  the  dipolar  geomagnetic  system  for 
use  with  the  ray  trace  program.  The  location  of  the  trough  wall 
is  assumed  to  be  composed  of  the  great  circle  segments  which 
connect  the  input  points.  Variations  in  the  plasma  frequency 
in  the  horizontal  direction  then  become  a function  of  the  distance 
from  the  trough  wall.  Distance  is  taken  to  be  positive  north 
of  the  wall  base  and  negative  for  points  to  the  south  The 
trough  wall  used  for  the  current  study  is  shown  in  Figure  1. 

Details  of  the  computer  program  which  produces  the  coordinates 
for  the  trough  wall,  program  WALL,  are  given  in  Appendix  A. 

The  parameters  which  vary  with  the  distance  from  the  trough 
wall  are  the  F-layer  and  E-layer  maximum  electron  density  values 
These  values  are  given  in  terms  of  critical  frequencies  (foE,foFt;  = 
maximum  plasma  frequency  of  the  respective  layer)  or  in  N^max 
(e lec t rons /cm  ' for  the  respective  layer  maximum),  where  the  two 
terms  are  related  by  the  formula; 

N max  = 12‘f00 . • fo“^ 
e 

The  foE  value  varies  from  0.1  MHz  south  of  the  trough  wall  to  2,5 
MHz  north  of  the  trough  wall;  the  height  of  the  maximum  electron 
density  (h  max)  is  assumed  to  be  120  km.  The  foF2  values  north  and 
south  of  the  trough  wall  are  input  quantities  to  the  ray  trace 
program  as  is  the  width  (or  steepness)  of  the  wall. 

Variation  across  the  trough  wall  is  taken  to  be  linear  in 
the  square  ol  the  plasma  frequency  which  is  equivalent  to  a linear 
variation  in  the  electron  density.  The  height  of  the  F,/  layer 
maximum  is  taken  to  be  constant  at  y50  km.  Electron  density  variation 
across  the  trough  wall  is  given  in  Figure  2 
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Figure  . . Variation  ol  loE  and  loF^  across  the  Trough  Wall 


Here  is  the  foF.  (in  MH2 ) at  the  bottom  of  the  trough 

and  foF2  is  the  foF^  at  the  top  of  the  trough.  W is  the 
width  of  the  trough  wall  in  degrees. 

The  various  trough  wail  mode  is  examined  in  this  study  are 
listed  in  Table  1.  They  were  derived  from  a trough  wall  model 
determined  by  Pike  ( 1 'Yu ) from  Alouette  topside  ionograms . 


Table  1.  Trough  Wall  Models  Used 

foF,.  foF.  . 

max  min  W 


Model  I is  the  median  model,  II  and  III  are  derived  from 
the  upper  and  lower  quartiles  of  the  model  gradient,  while 
models  IV  (most  shallow  trough),  and  V (steepest  and  widest  trough 
wall)  are  extremes  determined  from  the  model. 


The  vertical  variation  of  the  plasma  frequency  is  broken 
into  two  segments.  Up  to  the  E-layer  maximum  at  IkO  km.,  it  is 
composed  of  a Chapman  layer  with  a scale  height  of  10  km.  From 
120  km.  to  350  km.  the  variation  is  sine  squared  in  the  square  of 
the  plasma  frequency,  i.e. 
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where 


o _ II  *^maxE  _ 1? 

2 h ^-h  ^ " 2 

maxF  maxE 

h is  the  height  in  km. 


h-12o. 

230. 


Sample  vertical  profiles  are  given  in  Figure  . Details  on 
the  electron  density  model,  subroutine  BCHAP,  are  given  in  Appendix  B. 


The  sine  squared  segment  from  120  to  j50  km.  was  used  because 
it  was  felt  that  this  is  a more  realistic  presentation  of  the 
polar  ionosphere.  This  approach  made  the  model  less  sensitive  to 
variation  in  the  height  of  the  F2  layer  than  it  would  be  in  a 
Chapman  model.  Detailed  comparisons  of  the  two  models  are  given 
in  Section  3*£-  A sketch  of  the  two  types  of  profiles  is  given  in 
Figure  4. 


Figure  ••  Comparison  of  Sine  Squared  and  Chapman  Vertical  Profiles 
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. Displays  Used 

Since  the  quantity  which  the  OTH  backscatter  system  yields 
as  an  indication  of  target  range  is  the  group  delay,  the  group 
delay  must  then  be  converted  to  ground  range.  As  a possible 
model  for  determining  ground  range  from  group  path  length,  the 
sine  squared  model  with  an  F-layer  maximum  of  3.O  MHz  at  km 

and  an  E-layer  maximum  of  01  MHz  at  120  km.  was  used  The  curve 
of  group  path  versus  ground  range  for  a transmission  frequency 
of  5 MHz  is  given  in  Figure  5- 

To  get  a measure  of  the  deviation  in  azimuth  angle  and  error 

in  ground  range  caused  by  the  trough  wall  model,  ray  tracing 

results  from  the  disturbed  models  previously  described  were  compared 

with  results  from  the  horizontally  stratified  MHz  model  described 

above,  (j  MHz  was  considered  a reasonable  choice,  since  it  is 

well  within  the  range  of  foF2  values  to  the  north  of  the  pole- 

ward  wall  model,  as  shown  under  foF2  in  Table  1.  It  also 

max 

is  a reasonable  value,  if  the  standard  ITS-yft  foF2  map  would  be 
considered  the  basis  for  the  computation  of  radar  range.)  The 
comparison  of  range  and  azimuth  results  for  computations  using 
the  realistic  and  the  horizontally  stratified  models  was  made 
for  several  frequencies  (i.  to  9 MHz),  a selected  azimuth  range 
of  transmissions  (20.^*^  to  ■•2^  East  of  North  in  the  geographic 
coordinate  system),  and  for  elevation  angle  from  0°  to  a maximum 
o f . 

The  procedure  for  computing  range  errors  was  to  consider  as  a 
reference  the  range  versus  group  path  derived  from  the  j MHz  model 
as  given  in  Figure  instance,  for  a transmission  frequency  of 

MHz,  an  elevation  angle  of  and  an  azimuth  angle  of  the  model 

referred  to  as  Trough  I gave  a group  path  length  of  km.  and  a ground 

range  of  229^'  km.  The  MHz  standard  model  would  yield  a ground  range 

of  2^D‘*  km.  for  a group  path  length  of  2ohh  km.  meaning  that  the  ground 
range  of  the  target  determined  from  a horizontally  stratified  model  was 
actually  269  km.  longer  than  the  actual  ground  range.  Thus  the  range 
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error,  which  would  be  observed  if  a realistic  model  were  not 
available  and  radar  range  computations  were  based  on  a (reasonable) 

3 MHz  stratified  model,  would  be  +2fa9  km. 

A similar  procedure  was  followed  in  estimating  the  power 
difference  for  the  return  from  a target  whose  path  was  influenced 
by  the  trough  against  that  of  a target  that  was  observed  via  a 

MHz  horizontally  stratified  ionosphere  (outside  the  trough  wall), 
using  the  3 MHz  standard  ionosphere.  The  power  level  loss  was  computed 
as  a function  of  group  path  length  as  shown  in  Figure  6.  The  power 
level  for  any  given  trough  model,  frequency,  elevation,  and  azimuth 
angle  was  then  compared  against  the  standard  according  to  its  group 
path  length.  In  this  way  it  was  determined  whether  the  target 
would  appear  stronger  or  weaker  as  a result  of  trough  bending,  and 
by  how  much  Details  on  the  power  calculations  are  given  in 
Section  3. 2. 

j 1 Comparison  of  Sine  Squared  Versus  Chapman  Profile 

Before  using  the  3.O  MHz  sine  squared  layer  with  its 
maximum  at  39O  km.  as  the  standard  for  determining  the  ground 
range  and  power  as  a function  of  group  path,  the  curves  were 
also  computed  for  foF2  = 2.5  and  3.5  MHz  (Figures  7 and  lO)  and 
for  the  . MHz  case  with  the  maximum  at  3OO  and  1-00  km.  (Figures  8 
and  11). 

When  the  value  of  foF2  was  varied,  the  skip  distance  varied 
considerably  as  could  be  expected,  but  the  ground  range  derived  from 
a given  group  path  varied  by  less  than  50  km  with  the  change  of  foF2 
from  . .5  to  3.5  MHz.  The  power  values  for  a given  group  path  length 
varied  by  5 dB  for  values  away  from  the  skip  distance. 

When  the  height  of  the  F layer  maximum  was  varied  for  the  case 
of  • . ) MHz,  the  skip  distance  varied  little.  The  ground  range  for  a 
given  group  oath  varied  by  less  than  3O  km  in  most  areas  and  the 
power  loss  varied  by  less  than  2 dB  for  elevation  angles  below  the 
skip  elevation  angles. 
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TRANSMISSION  FREQUENCY  IS  5 MHz. 


Power  Loss  In  Decibels 


3.0  MHz 


Group  Path  Length  in  km 

FIGURE  10.  POWER  LOSS  VERSUS  GROUP  PATH  LENGTH  FOR  SINE  SQUARED 
PROFILES  WITH  foF2  = 2.5,  3.0,  AND  3.5  MHz  LOCATED  AT 
350  km  ALTITUDE 


>^^350  km 

f 400  km 
300  km 


Group  Path  Length  in  km 


FIGURE  11.  POWER  LOSS  VERSUS  GROUP  PATH  LENGTH  FOR  SINE  SQUARED 
PROFILES  WITH  foF2  = 3.0  MHz  LOCATED  AT  300,  350,  AND  400  km 


To conputv L he  results  derived  from  the  sine  squared  model  with  those 
trom  a Chapman  layer  model,  rays  were  traced  through  a j.O  MHz  Chaoman 
layer  with  scale  height  of  o7.^  km  and  with  the  maximum  at  heights  of 
, i‘^0 , and  400  km. 

The  variation  of  ground  range  for  a given  group  path  was  in  the 
order  of  100  km  over  the  three  models  (Figure  9)^  three  times  as 
sensitive  to  a 100  km  change  in  h as  the  sine  squared  layer. 

The  case  for  the  maximum  at  jOO  km  closely  followed  the  j .0  MHz  sine 
squared  standard  which  was  used.  The  power  loss  variation  for  a 
given  groun  path  varied  by  as  much  as  8 dB  and  was  considerably 
ditferent  in  apnearance  from  the  sine  squared  curves  (Figure  12). 

Again  the  km  case  most  strongly  resembled  the  sine  squared 

s t andard  . 


The  chosen  sine  squared  model  thus  produces  results  which  are 

less  sensitive  to  h or  foF^  changes  than  those  derived  from  a 
max 

Chanman  model. 

..  Power  Loss  Calculations 

Power  loss  calculations  were  for  backscattered  power  returned 
from  ground  reflections  lor  the  area  covered  by  a beam  l'^  in 

, , 1 o . . , 

elevation  and  1 in  azimuth. 


The  power  returned,  P^ , is  given  by: 


' r i-TT 


'■V  ‘T 


where  P is  the  transmitted  power, 
o 

The  quantity  used  as  the  power  loss  in  db  is  actually 


U log^^  (^)  = 1 + 2U  log^  ^ + 1C  log^JJ^+  10  log^^C-^) 


o 

p 

' A 

is 

tho 

wave  length  in  km. 

A 

is 

the 

ray  spreading 

P 

0~ 

is 

the 

backscatter  cross-section 

p 

w 

is 

the 

pulse  width  (Ac  km.  was  used) 

is 

the 

rms  group  path  deviation. 

I 


Power  Loss  in  Decibels 


FIGURE  12.  POWER  LOSS  VERSUS  GROUP  PATH  LENGTH  FOR  CHAPMAN  PROFILES 
WITH  foF2  = 3.0  MHz  LOCATED  AT  300,  350,  AND  400  km 


1- 


Other  power  loss  factors  which  are  not  included  at  the  present 
time  are:  absorption,  transmitter  and  receiver  gain.  Absorption 
was  not  considered  in  these  model  computations,  since  it  would  not 
impact  differently  on  rays  propagated  via  a disturbed  or  a smooth 
ionosphere.  The  transmitter  and  receiver  antenna  patterns  may  be 

readily  imposed  as  a function  of  frequency,  azimuth,  and  elevation 
on  the  results  reported  in  Section  U.  A more  detailed  discussion 
of  the  computation  of  backscatter  and  spread  losses  may  be  found 
in  Chapter  IV  of  RADC-TR-77-6O . 
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‘-t  .0  Signiticaiit  Observations 


Rays  wore  traced  through  the  trough  models  given  in  Table  1 
and  the  power  ditferences,  range  errors,  and  azimuth  deviations 
were  computed.  Transmission  frequencies  used  were  hj  o,  'f , 
and  1 MHz.  It  was  found,  for  the  models  used,  that  rays 
escaped  for  most  elevation  angles  for  frequencies  greater  than 
t MHz.  Trougti  IV.  the  model  with  the  weakest  gradient  and  lowest 
foFv  poleward  of  the  trough  wall,  also  gave  poor  results  for  fre- 
quencies of  5 iitid  o MHz.  Rays  were  traced  every  degree  in 
elevation  from  0 to  ^0^  and  for  azimutli  angles  from  c'u . ° to 

in  the  geographic  coordinate  system.  The  hypothetical  trans- 
mitter was  located  at  -.a.h  N geographic  latitude  and  o7  • '*^W 
geographic  longitude.  The  dipolar  geomagnetic  pole  was  assumed 
located  at  c N latitude  and  ■.  '.7'-l  W longitude.  Azimuth 
angles  of  -1^  and  greater  were  unaffected  by  tlie  trough  wall. 

Using  Trough  Models  II  (weakest  gradient,  lower  quartile),  1 
(median)  and  111  (strongest  gradient,  upper  quartile),  the  effects 
of  an  increase  in  the  trougli  wall  steepness  on  focusing  and  defocusing 
can  be  shown.  Figure  1 ■ a,  b,  and  c shows  the  difference  of  the 
backscatter  power  derived  from  the  realistic  trough  wall  model  and 
that  derived  from  the  horizontally  stratified  ■ MHz  model  as  a 
function  of  elevation  angle  and  azimuth  angle  of  the  launched  ray. 

In  general,  the  power  compared  to  the  - MHz  standard  ionosphere  is 
enhanced  over  a wide  region  for  Trough  II  while  this  enhanced  region 
becomes  increasingly  narrow  in  azimuth  as  the  trougli  wall  steepens  to 
Trough  III  and  an  ever  increasing  area  of  power  degradation  appears 
although  most  relative  losses  are  less  than  10  dB. 

Relative  enhancement  of  power  is  greatest  for  those  azimuth  angles 
where  the  rays  just  graze  the  trough  wall.  This  is  in  part  due  to 
the  fact  that  rays  heading  more  directly  into  the  trough  wall  are  bent 
down  fairly  sharply  giveing  a group  path  length  which  is  nearer  to  the 
■ .0  MHz  skin  distance.  Even  though  these  backscatter  signals  do  have 
a power  comparable  to  the  more  grazing  rays,  they  show  less  power 
returned  if  compared  with  the  horizontally  stratified  model. 
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th  in  Degrees 


TROUGH  AND  GRADIENT  RANGE 


I 


Rays  which  escape  the  eltects  of  the  trough  wail  (those 
having  azimuth  angles  greater  than  approximately  -l^)  experience 
an  apparent  enhancement  of  power  if  the  j;,MHz  reference  model  is 
used.  This  occurs  because,  in  fact,  these  rays  travel  through 
a d MHz  horizontally  stratified  ionosphere.  In  order  to  eliminate 
this  false  enhancement  a 2 MHz  reference  is  used  for  these  rays 
and  for  rays  which  may  be  slightly  influenced  by  the  trough  wall 
but  which  for  the  most  part  travel  in  the  2MHz  Ionosphere.  Thus 
for  all  rays  with  azimuth  angles  greater  than  the  2 MHz 

reference  is  used. 

The  range  and  azimuth  errors  shown  in  Figures  1';  and  16 
show  a behavior  similar  to  one  another  in  that  large  range  errors 
generally  occur  for  the  same  azimuth  and  elevation  angles  as  large 
azimuth  deviations.  Range  errors  decrease  with  increasing  wall 
steepness  mainly  due  to  the  fact  that  rays  are  bent  more  sharply 
downwards  and  hence  travel  less  far  in  a direction  tangential  to 
the  transmitter.  The  azimuthal  deviation  also  decreases  with  in- 
creasing wall  steepness  uptoapoint  although  rays  are  bent  sideways 
more  sharply  as  the  steepness  increases.  This  is  again  due  to 
the  fact  that  rays  travel  farther  after  they  are  bent  by  a weak  trough 
wall.  However,  as  the  trough  wall  continues  to  increase  in  steepness, 
although  rays  travel  less  far  after  bending,  they  are  bent  so  sharply 
sideways  that  the  overall  azimuth  deviation  increases. 

This  behavior  can  be  more  clearly  understood  by  looking  at 
the  ray  paths  as  a function  of  range  and  azimuth  from  the  trans- 
mitter which  are  shown  for  Troughs  II,  I,  and  III  in  Figures  17, 

18,  and  19  respectively.  The  elevation  angle  in  these  figures 
is  7° • For  the  azimuth  angle  the  azimuth  deviations,  range 

errors,  and  total  ground  range  are  given  in  Table  2. 
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Elevation  in  Degrees  Elevation  in  Degrees  Elevation  in  Degrees 


b)  Trough  I 


azimuth  in  degrees 


c)  Trough  III 


azimuth  in  degrees 


FIGURE  15.  AZIMUTHAL  DEVIATION  AS  A FUNCTION  OF  ELEVATION 
AND  AZIMUTH  FOR  5 MHz 


Elevation  Angle  in  Degrees  Elevation  Angle  in  Degrees  Elevation  Angle  in  Degrees 


v<%%vi:’v9 


m 


IMl 


Range  from  the  Transmitter  in  km 


Plasma  Frequency  Contours  (MHz) 
Are  For  Height  = 350  km 


3000  t III 


Transmi tter 

FIGURE  18.  RAY  PATHS  PLOTTED  AS  A FUNCTION  OF  RANGE  AND  AZIMUTH  FROM 
THE  TRANSMITTER  FOR  7"  ELEVATION  ANGLE,  TRANSMISSION 
FREQUENCY  5 MHz,  AND  IONOSPHERIC  MODEL  TROUGH  1 


Range  from  the  Transmitter  in  km 


Plasma  Frequency  Contours  (MHz) 
Are  For  Height  = 350  km 


Transmi tter 


FIGURE  19.  RAY  PATHS  PLOTTED  AS  A FUNCTIOf!  OF  RANGE  AND  AZIMUTH  FROM 
THE  TRANSMITTER  FOR  7'’  ELEVATION  ANGLE,  TRANSMISSION 
FREQUENCY  6 MHz,  AND  IONOSPHERIC  MODEL  TROUGH  III 


Table  2. 


Comparison  of  Azimuth  Deviations  and 
for  Three  Trough  Models  at  Elevation 
Azimuth  Angle  - and  Frequency  - 


Ra[ige  Errors 
Angle  - 7^, 
MHz. 


Az  imu  t h 
Deviation  in 
Mode  1 degrees 

II  -5.^ 

1 - .o 

HI  -5.  < 


Range 

Ground 

Error  in 

Range 

km . 

km . 

(500 

23--' 

270 

-115 

270 

1 'o2 

Ground  Reflection 
Az imuthal 
Angle 

50.88 

51.-7 

53.‘-V 


The  ground  reflection  azimuthal  angle  given  in  Table  e is  the 
azimuth  angle  East  of  North  geomagnetic  at  which  the  ray  was 
traveling  when  a ground  reflection  occurred  In  this  case, 
while  the  ray  for  Model  II  was  bent  to  an  azimuth  angle  less  than 
Model  I,  its  azimuth  deviation  was  greater  due  to  the  increased 
ground  range  gained  (particularly  that  gained  after  off  great 
circle  bending).  However,  for  Model  III  the  azimuth  angle  to 
which  the  ray  was  bent  was  great  enough  that  there  was  an  increase 
in  azimuth  deviation  over  Models  II  and  I even  though  the  ray 
went  less  tar. 


.0  Conclusions 


For  the  great  majority  ot  the  rays  examined  the  power  return 
t rom  targets  located  at  their  landing  points  should  be  great  enough 
to  allow  detection,  it  detection  at  the  same  trequency  under  smooth 
ionospheric  conditions  would  have  been  possible.  The  assumptions 
made  about  the  target's  position,  it  a trough  model  is  not  considered 
in  such  a determination,  could  be  ott  in  range  by  up  to  j>0U  km  and  !:ave 
an  azimuth  deviation  ot  up  to  (the  equivalent  ot  cL'JO  km  at  the 
ranges  we  are  considering)  It  a trough  model  is  considered,  target 
positions  are  still  strongly  dependent  on  the  chosen  trough  wall 
gradient.  Ranges  vary  trom  +1U0  to  - )0  km  relative  to  ranges 
derived  trom  the  median  model,  it  the  gradient  is  changed  to  the 
lower  quart  ile  or  the  upper  quart  ile  value  respectively. 

There  is  also  a region  close  to  the  trough  wall  boundar>  where 
multi-path  occurs  thus  giving  rise  to  tlie  appearance  ot  two  targets 
where  one  actually  exists. 
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APPENDIX  A 


PML  152 


NAME:  WALL,  revision  0,  program,  PML  152 

CATEGORY:  Preprocessing  program  for  use  with  ray-tracing  programs 

TITLE:  Trough  wall  preprocessing  program  for  subroutine  BCHAP 

LANGUAGE:  CDC  extended  Fortran  - version  h 

PROGRAMMER:  B.M.  Langworthy,  Parke  Mathematical  Laboratories,  Inc. 

DATE:  October  21,  I976 

DESCRIPTION 

Program  WALL  accepts  coordinates  of  a path  in  the  geographical 
coordinate  system  which  specifies  the  poleward  side  of  the  arctic  trough 
bottom.  It  then  converts  these  points  to  the  dipolar  geomagnetic  system 
and  computes  certain  quantities  needed  in  the  computation  of  plasma 
frequency  by  ray-trace  subroutine  BCHAP.  These  quantities  are  put  out 
on  a tape  or  permanent  file  for  use  by  the  ray-trace  program. 

INSTRUCTION  SET 

To  use  program  WALL,  the  following  data  must  be  entered: 

Card  1 (AIO) 

NAME  - This  is  a 10-character  alphanumeric  name  for  the  trough 

wall.  It  will  be  used  by  subroutine  BCHAP  as  an  electron 
density  subroutine  name  and  will  appear  on  each  page  of 
ray-trace  output.  This  name  must  be  punched  in  columns 
1 through  10. 

Card  2 (2F10.2)  GLAT(1),  GL0N(1) 

GLAT(1)  is  the  geographic  latitude  in  degrees  of  the  first  point. 

GL0N(1)  Is  the  ge  -graphic  longitude  in  degrees  of  the  first  point. 

GLAT  must  be  punched  in  columns  1 to  10  and  must  contain  a 
decimal  point.  GLON  must  be  punched  in  columns  11  to  20  and 
must  also  contain  a decimal  point. 

Carcfe  and  on 

contain  subsequent  points  along  the  trough  wall  and  data  should 
be  entered  in  the  same  manner  as  on  Card  2. 

.^(Eur.DI  0 ;->..DE  PL-IIK-.XT  FILI.iD 
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Output  from  program  WALL  includes  printed  output  and  an  unformatted 
tape  or  permanent  file  refered  to  as  TAPE6.  TAPEd  is  described  under 
FILE  DESCRIPTIONS.  A sample  of  the  printed  output  is  shown  in  Attachment  1. 
For  each  input  point  the  following  information  is  given: 

GLAT  - the  geographic  latitude  of  the  point  in  degrees 
GLON  - the  geographic  longitude  of  the  point  in  degrees 

(?  - the  dipolar  geomagnetic  colatitude  in  degrees 

- the  dipolar  geomagnetic  longitude  in  degrees 
a - the  length  in  degrees  of  the  great  circle  segment  connecting 
( , (f*  ) with  the  next  point 

3 - the  azimuth  angle  in  degrees  of  the  great  circle  segment 
connecting  ( ^ , ) with  the  next  point 

(f  - one  half  the  angle  in  degrees  formed  by  the  two  great  circles 
which  meet  at 

The  message  "TAPEo  HAS  BEEN  WRITTEN"  will  appear  after  the  table  to 
indicate  that  ^he  program  has  completed  writing  TAPEb. 


STORAGE  REQUIRED 

Program  WALL  requires  40000  octal  words  of  core  storage. 
ALGORITHM 


Points  inputted  in  the  geographic  coordinate  system  are  converted  to 


points  ( O where: 


, Lh 


entrv 


(5  ^ is  the  dipolar  geomagnetic  colatitude  of  the  i 
c|'  ^ is  the  dipolar  geomagnetic  longitude  of  the  i^*^  entry 


Several  other  quantities  are  computed  at 
electron  density  subroutine  BCHAP.  They 
defined  and  illustrated  on  the  following 


this  point  which  are  useful  for 

are  a . , ^ , and  6 i which  are 
1 ■ 1 

page  . 


L. 


'1 
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North  pole 

« 


-A 


a = cos  ^(cos  G.cos  (?.  + sin6>.sin6».  ,cos(cf  - (T  )) 

1 '11+1  i i+1  ^ Vi+1  ^ i' ' 


^ , cos  O , , - cos  (•■'.cos  a 

li  . -1  / 1+1  1 i. 

1 = cos  ( : = i-J 


If  i > 1 


s in  (h  , s in  a . 

i 


.i  = cos  ^(cosQ  . ,cosG.  1 +sint>.  .sin  6^.  ,cos(;:.''  )') 

^•-<1  J--I  1+1  1-1  1+1  ' > 1+1  ^ i-lv 


, cos  a . - cos  a,  .cos  a 

ri  = 1/2  cos'^  ( i) 


^ . = cos 

• lit 


sin  a,  ,sin  a. 

1-1  1 

cosG  , cos  a. 


^,i\ 

\ sin^,  ^_j^sin  a^^  / 

If  i = 1, 

"If  i = NENT  or  the  last  entry,  = 1 radian,  ~ ^ i = ■—  . 


SPECIAL  CAITK^JS  AND  FEATURES 


The  number  of  point  entries  should  not  e.xceed  oO. 


TLMING 

Execution  time  is  .05  seconds  for  10  points 


error  messages 


THERE  IS  NO  INPUT  DATA. 

This  message  will  be  printed  in  data  cards  if  point  entries  are  inadvertently 
omitted . 

SUBROUTINES 

DICOORD  - The  subroutine  converts  points  from  the  geographic  to  the 
dipolar  geomagnetic  system. 

ACCURACY 


The  accuracy  of  the  trough  wall  outputted  by  WALL  is  dependent  on 
the  accuracy  and  closeness  of  the  input  points. 


FILE  DESCRIPTIONS  WALL( INPUT, OUTPUT, TAPEL =INPUT, TAPEo ) 

Program  WALL  writes  data  onto  TAPe6  to  describe  the  trough  wall  to 
Subroutine  BCHAP.  TAPEo  contains  the  following  data  in  unformatted  form. 


Record  1 


Record  2 


Record  j 
Record  4 


Record  5 
Record  6 


NAME,  NE  NT 

NAME  - an  alphanumeric  rame  for  the  wall  data  set 
NENT  - the  number  of  entry  points  (<  30) 

(DCLAT(i),  i = 1,NENT) 

DCLAT  - array  of  dipolar  geomagnetic  colatitudes  in  radians 
(DLON(i),  i = 1,NENT) 

DLON  - array  of  dipolar  geomagnetic  longitudes  in  radians 
(SDA(i),  i = 1,NENT) 

SDA  - length  of  great  circle  segments  connecting  points 
in  radians 

(BETA(i),  i = IjNENT) 

BETA  - the  azimutli  angle  of  the  line  segment  in  radians 
(CAM(i),  i = 1,NENT) 

CAM  - one  half  the  angle  between  consecutive  segments 


in  radians 
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DECK  SETUP 

LAiWG,T10,CM50000.  NO.  NAME 

FTN. 

REQUEST(TAPEo,*PF) 

LGO. 

GATALOG ( TAPEo , UTODWALLX369381 6, ID  =LANGW , RP  =999 ) 

7/8/9 

Deck  of  program  WALL  and  subroutine  DICOORD 

7/8/9 

Data  cards 

d/7/6/9 
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t 

t 

I 

t 

t 

[ TPOU.ri  WALL  C6JT 


GEOG 

ULOG 

L- 1 r- 

I'l  7 

L AT  I 1 Ub£ 

L ONG. 

cn.Aj . 

1 OMG. 

*.)  j£  A 

G A r '1 A 

51. 50  0 

?95. og: 

0 - . ^ 1 3 

c . 0 ^ 

20' 

7i. 11. 

9 . 
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w W V • . . 

1-  c . 1.  3 L. 

io.‘.7- 

0.  ?3" 

! C • V 1 7 

t .12" 

53.  70  C 

3G5. 000 

cl. , 3 fc  0 

2 ; . 5 9 u 

3 . 25  , 

^ f _ 

c ’ .1 67 

5 5 , H ij  C 

3 10 . ^ , 

C ^ c 

.6.: 1; 

i . Ju  i 

t *-  • L 3 

5 , . t b 

57.  ^0  0 

315. 0G5 

Zc,  1.  c 7 

3i  .71  i 
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7 • 7 3 > 

*■  ’.135 

55.  LG  V, 

^ ^ M • V ' 

cl.. 7-^ 

. 1 5 c 

3.  35*- 

t.  L . :7^ 

t.  0 .-.5- 

D 1. 90  0 

f pc  r, 

<->  C ^ 9 J J - 

1 i. 723 

0 ?.??:’ 

3. 

t.  J • L L- 

5 •:  . •”  8 ’ 

0 9 . 3 J c 

3 w • U C 

1 0 . 0 c 3 

1 1 . 1 3 _ 

2 . j s . 

- • 7 r A 

5 ■'.221 

0 tf . 30  w 

3 35.  0 00 

1 - 3 ^ 

7 J .651 

2.70- 

7^. 

c ’ . 33, 

6 S . (^0  L 

W H w • c ^ ■ 

. c . / 1 

7 -y  . 7 6 1 
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^ • V w 

-j  • i L . 

TAPtb  HAS  3E£N  WPITTE:.. 


Attachment  1.  Sample  Output  from  Program  WALL 
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PMLlSl 


NAME : 

CATEGORY : 

TITLE: 

LANGUAGE: 

PROGRAMMER: 

DATE: 


BCllAP.  revision  0,  subroutine,  PML  l^i 
Electron  density  subroutine  for  Ray-Tracing  Program 
Double  Chapman  Layer  for  Trough  Simulation 
CDC  Extended  Fortran  - version  h 

B.M.  Langworthy,  Parke  Mathematical  Laboratories,  Inc. 
October  21,  l97b 


DESCRIPTION 


The  electron  density  model,  BCllAP,  is  a replacement  for  electron  density 
subroutine  CHPTRH,  PML  L2&.  Whereas  in  subroutine  CHPTRH,  the  trough  wall 
was  located  along  a line  of  constant  dipolar  geomagnetic  colatitude,  subroutine 
BCllAP  can  use  an  arbitrary  trough  wall  orientation.  The  electron  density 
model  is  a double  Chapman  layer  with  variations  as  a function  of  perpendicular 
distance  from  the  trough  wall.  At  present  the  only  Chapman  layer  parameters 
which  vary  are  the  cr it ical frequonc ies  of  the  two  layers.  The  layer  height 
and  thickness  of  both  layers  remain  fixed.  The  location  of  the  trough  wall 
and  other  parameters  needed  by  subroutine  BCllAP  are  supplied  on  input  TAPE6. 
This  file  is  created  by  program  WALL,  PML  1‘2. 

INSTRUCTION  SET 

The  use  of  CHPTRH  is  similar  to  that  of  other  electron  density  subroutines 
used  by  the  ray-trace  program.  In  the  event  that  the  selective  load  method 
described  in  PML  121  and  AFCRL-TR-75"Ojl9  is  being  used,  the  SLOAD  card  will  be: 

SL0AD(BSUBS,  BCIL^  P,  and  any  desired  magnetic  field,  collision 
frequency,  and  perturbation  subroutines) 

No  auxiliary  subroutines  are  required  by  subroutine  BCllAP. 

Input  required  by  the  subroutine  is  in  two  forms.  First,  Che  location 
of  base  of  the  trough  wall  is  defined  by  coordinates  _,iven  on  TAPED.  This 
information  is  explained  in  detail  under  FILE  DESCRIPTIONS.  For  further 
infonnation  on  it,  refer  to  PML  1A2  which  describes  the  process  of  generating 
TAPE6.  The  other  information  required  by  subroutine  BCllAP  is  input  to  the 
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W-array  ot  the  ray-trace  program.  Three  locations  must  be  defined  as 
tOi lows : 


W-array 
locat ion 

101 

102 


Quant ity 

F2MX  The  maximum  plasma  frequency  of  the  F2  layer 

in  MHz  at  the  top  of  the  trough  wall. 

F2MN  The  maximum  plasma  frequency  of  the  F2  layer 

in  MHz  at  the  bottom  of  the  trough. 


lOj  PMAX  The  width  in  radians  of  the  trough  wall.  If 

PMAX  is  entered  in  degrees,  a "1"  must  be  placed 
in  column  16  of  the  input  card. 


Subrout ..ne  BCHAP  uses 


the 


following  infonr.ation  from  the  ray-trace  program: 


R(l)  - r - 
R(2)  - O - 
R(3)  - 0 - 
E ART  HR 
F 

PI  - rr  - ] 
PIT2  -2  ii-y 

PID2  - ij  -j 

2 


radius  in  km. 

dipolar  geomagnetic  colatitude  in  radians, 
dipolar  geomagnetic  longitude  in  radians. 

(also  W(2))  the  radius  of  the  earth  in  km. 
(also  W(3))  the  transmission  frequency  in  MHz. 

constants . 


The  subroutine  then  returns  the  following  quantities: 


MOUX(l) 

X 

PXPR 

PXPTH 

PXPPH 

HM/\X 


- an  alphanumeric  identifier  which  is  supplied  by  TAPEb. 
identifying  the  location  of  the  trough  wall. 

- refractivity  due  to  plasma  frequency  (dimensionless). 

- partial  derivative  of  refractivity  with  respect  to  the 
radius  (/km). 

- partial  derivative  of  refractivity  with  respect  to 
colatitude  (/radian). 

- partial  derivative  of  refractivity  with  respect  to 
longitude  (/radian) 

- the  maximum  iieight  of  the  electron  density  in  km. 

(This  is  used  to  determine  when  penetration  has  occured.) 
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STOR/V-c  REQUIRED 

Subroutine  BCHAP  requires  7OO  octal  words.  When  used  in  the  ray-tracing 
program  with  index  of  refraction  AVWFNC  and  magnetic  field  DIPOLY,  5OOOO  octal 
words  are  required. 

ALGOP  ITI{M 


The  plasma  frequency  (or  electron  electron  density)  model  given  by 
BC51AP  can  be  broken  into  three  separate  computational  areas; 

1)  computation  of  wall  location  parameters;  2)  computation  of  perpendicular 
distance  from  the  wall  and  angular  components  relative  to  the  wall;  and 
y)  computation  of  the  trough  wall  and  the  vertical  profile.  The  first  area 
of  computation  is  performed  by  program  WALL  but  will  be  briefly  presented 
here  since  the  definition  of  terms  will  be  needed  subsequently. 


1)  Computation  of  wall  location  parameters 


The  great  circle  segments  which  give  the  location  of  the  trough  wall  are 
specified  by  their  endpoints  in  the  dipolar  geomagnetic  system.  The  i th 
endpoint  will  be  referred  to  by  wheret?is  the  dipolar  magnetic 

colatitude  and  j/'  is  the  dipolar  magnetic  longitude.  The  point  ) 

without  subscripts  is  the  point  for  which  the  plasma  frequency  and  its 
derivatives  are  to  be  computed.  The  geometry  of  the  situation  is  given 
below.  ' 


\ 


The  quantities  which  will  be  needed  for  later  computations  are: 


a . 
1 


the  gri_at  circle  path  length  in  radians  between  points  i and  i + 1, 
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the  clockwise  positive  angle  in  radians  with  vertex  at  point  i 
between  due  north  and  the  path  connecting  points  i and  i+1 . 
one  half  the  angle  in  radians  formed  by  points  i-1,  i,  and  i+1. 

Excluding  exceptions  for  end  points, 


a.  =cos  ^(cost^.cost^,  , +sint5,sinc^.  ,cos(d'.  , ~ ^ ■)) 
1 ' 1 1+1  1 1+1  ^ I 1+1  11'' 


- cos  O . cos  a . 
1+1  1 1 


s in  <L'^  . s in  a , 
1 1 


cos  a - cos  a.  .cos  a. 

r.  = 1'2  cos-'/ i 

' s in  a . , s in  a . 

1-1  1 


where  cos  a^  =cos(9,  ,coS(f>.  , +sin^.  ,sini3,  ,cosl  Q\  - 4' . 

1-1  1+1  1-1  1+1  ^'1+1  f i-l' 

2)  Computation  of  the  perpendicular  distance  from  the  wall  to  point  {&,(]>) 

To  eliminate  any  ambiguity,  perpendicular  distance  from  the  wall,  P,  will 
be  considered  positive  if  the  point  ( ff , (^  ) is  on  the  north  poleward 
side  of  the  wall  and  negative  if  it  is  on  the  other  side.  Since  all 
points  except  those  lying  on  a given  great  circle  will  have  at  least 
two  perpendiculars  to  the  great  circle,  it  must  be  stated  that  the 
perpendicular  must  occur  with  the  segment  of  any  great  circle  which 
forms  a part  of  the  trough  wall.  There  may  also  be  cases  where  there  is 
no  such  perpendicular  due  to  "kinks”  in  the  wall  where  two  segments 
join.  The  first  problem,  then,  is  that  of  finding  the  proper  great 
circle  segment . 

To  find  P,  the  distance  from  the  curve,  and  7^,  the  angle  with  respect 

to  positive  ),  each  segment  is  tried  beginning  with  i = 1 

for  the  first  time  and  beginning  with  the  preceeding  i value  thereafter. 


- au  - 


f 


u'i*. 
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?>>  3-' 


19 


A, 

'fk 


t'i 


C - distance  from  (0,f)  to  {0-^,(p.) 

C = cos  (cosfV^cos6r  + s in  c?  ^s  iii^cos  ( (^  - .)) 

= angle  from  due  north  with  vertex  at  .)  to  (6~  , <j  ) 

y cosfv  - cos;>^cos  C 


J = cos  (- 


s in  ^>^s  in  C 


-) 


If 


< P ^ and  J-,  then  b = A - ' 1 

P is  positive  f Region  I 
proceed  to  * J 

If  p ^ i P i ~ - if  ^ *^*'“-*"  next  lowest  segment.^  Region  V 

< /O  ^ and  0 ^ - P '^  = ■’^i  ■ 

P is  negative  / Region  II 
proceed  to  * 1 


If 


If 


^ ^ ^ . and  y < fj  ^ ^ < ij:.  - 2-13' .,  then  B =./!?.  - ;3  ',1 

^ W.:.o  i 


P = c 

proceed  to  * 


Region  III 


It  ^ > pt,and  p - G - 2 J)  then  next  lowest  segmen^  Region  IV 


P =sin  ^(sin  B sin  C) 


J,  = distance  along  the  great  circle  connecting  ( C*  , (}>  ) and  (O  . ^ 

L ° ' 1’  I 1'  ' i+1’  T i+i 

to  the  perpendicular  from  ( £>■  , ) 


- 1 - 
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. -1 / cos  C » 

d.  = cos  1 1 

1 ' cos  r ' 

If  > a,  the  perpendicular  does  not  fall  within  the  desired  segment 

and  the  next  highest  (in  i)  segment  is  tried. 

If  the  perpendicular  occurs  below  the  point  the  segment  is 

taken  along  a constant  latitude  of  & Here  P “ <3  - O and 

, -1/  cos  C X 

d,  = - cos  ( ) 

■ ' cos  P ' 


To  find  , the  angle  between  the  perpendicular,  P,  and  the  positive  (9 
direction,  we  compute  two  additional  quantities. 

V : ^ ‘ 


■ 


{( <\) 


^ -1  j / 

y = cos  (-  cos(q)-('  .)  cos  <3  +Sin{{^-^  ^)s  in  cos  ^ 


It 


y = cos  ^(sin  B cos  d^) 

II  I 

V = - f ' - y 


For  the  region  below  y = fT 


- <*2  - 
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3)  Computation  of  the  trough  wall  and  the  vertical  profile 

The  trough  wall  currently  in  BCHAP  varies  as  a function  of  P,  the 
perpendicular  distance  on  inputted  curve,  in  the  maximum  values  of 
the  E-  and  F-  layer.  At  this  time  the  values  of  height  and  layer 
width  remain  fixed.  The  height  of  the  E-  layer  is  120  km  and  its 
scale  height  is  10  km.  The  height  of  the  F-  layer  is  350  km.  The  E- 
and  F-  layer  variations  are  given  below. 


f - plasma  frequency 
P 


Input  to  (F2MN)‘ 

W(102) 

.01MH“^ 


Input  to  w(lOl) 
(F2MX)^ 


6.25MH 


Y 


Input  to 
W(103) 


Given  P and  Y from  the  previous  computation  we  can  now  compute  the 
plasma  frequency  and  its  derivatives  at  the  point  (r,  O , (J'  ). 

h = r - r where  r is  the  radius  of  the  earth, 
e e 

If  70  < r < 120, 

f^  = f^  exp(l/2  (1  - Z,  - e ^1))  where  Z,  = 

pCj^  ' 1 ''  1 10. 


exp  (1/2(1-Zj^  - e 


•1 


'Hfb 


~ ? p 


s in  £>  s in  y 


If  120  < h < 350, 


where  H = — ,.~  120 . 


?(f?) 


JL  / .2  P 


? r 230  sin  H cos  H 

^(4)\ 

z'  p " -a  p ■"I  yT“  ■ ~JT~)  « 


a p ^ p 


^(<)  ^(C) 


"a  <=•>•  ? p 


* cos  "V 


7^  (f  ■ a>  p 


sin  <>  sin  V 


f = .01  if  P < 0 

Cf  — 


.01  + 6.24  if  0 < P < p 


0.2^  if  P > p 

— max 


= 6.24 


P P 

TT 


f^^  = (F2MN)^  if  P < 0 


= (F2MN)^  + (F2MX^  - F2MN^)  ^ of  0 < P < P 

P — n 

max 


= (F2MX)^  if  P > P 

— max 


^ ^^2  _ (r2>y^  - F2MN^l 

?)  p p 

''  • mav 


PML  151 


J 


PML  151 


aPECIAL  CAOTIONS  ANTD  FEATURES 

This  subroutine  as  all  other  electron  density  subroutines  for 
RAYTRACEBL  assumes  that  values  for  9 and  4 are  in  dipolar  geomagnetic 
or  computational  coordinates. 

If  the  number  of  entries  from  TAPEo  exceeds  30,  a new  dimension 
statement  must  be  inserted. 


TIM  I Nr: 


For  a large  number  of  rays  computed  using  BCHAP  execution  time 
was  about  1 second  per  ray.  This  will  aproach  1 .7  seconds  per  ray 
if  only  low  elevations  angles  are  considered.  Timing  was  done  on  1-hop 
rays  at  f 'equencies  from  5 to  9 MHZ. 


error  messages 


None. 


subroutines 

No  nonsystem  subroutines  are  required. 

accuracy 


The  accuracy  of  the  trough  in  BCIWP  is  dependent  on  the  closeness 
and  accuracy  of  the  points  describing  the  trough  wall  on  input  TAPEo. 


FILE  DESCRIPT IONS 


BCHAP  reads  data  off  of  TAPEo  during  program  initialization.  TAPEo 
contains  the  following  data  in  unformatted  form  generated  by  WALL(pmL  152). 
(Symbols  in  parentheses  relcr  to  symbols  used  in  the  algorithm  or.  pages  3-4 


1 


L 
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FILE  DESCRIPTIONS  (continued) 

Record  1 MODX(l),  NENT 

MODX(l)  - an  alphanumeric  name  for  the  wall  data  set 
NENT  - the  number  of  point  entries  (<  3O) 

Record  2 (XH(i),  i = 1,NENT) 

TH  - array  of  dipolar  geomagnetic  colatitudes  in  radians 

Record  3 (PH(i),  i = 1,NENT) 

PH  - array  of  dipolar  geomagnetic  longitudes  in  radians  ( ^'  . ) 

Record  4 (SDA(i),  i = 1,NENT) 

SDA  - length  of  great  circle  segments  connecting  points 
in  radians  ( a . ) 

Record  5 (BETA(i),  i = 1,NENT) 

BETA  - the  azimuth  angle  of  the  line  segment  in  radians  ( <^  . ) 

Record  d (GAM(i),  i = 1,NENT) 

GAM  - one  half  the  angle  between  consecutive  segments  in 
radians  ( 3^ 

DECK  SETUP 


Runs  using  BCHA.'  would  be  similar  to  other  ray-trace  runs  except 
that  a TAPE6  must  be  attached  and  the  SLOAD  card  given  on  page  1 should  he  used 


L. 


1 


